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Abstract 

Recent work on the spectrum of the Euclidean Dirac operator spectrum show that the exact 
microscopic spectral density can be computed in both random matrix theory, and directly from 
field theory. Exact relations to effective Lagrangians with additional quark species form the bridge 
between the two formulations. Taken together with explicit computations in the chGUE random 
matrix ensemble, a series of universality theorems are used to prove that the finite-volume QCD 
partition function coincides exactly with the universal double-microscopic limit of chUE random 
matrix partition functions. In the limit where Nf and iV c both go to infinity with the ratio Nf/N c 
fixed, the relevant effective Lagrangian undergoes a third order phase transition of Gross- Witten 
type. 
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1 Introduction 



Over the last five years it has gradually become clear that the eigenvalue spectrum of the Dirac 
operator in Euclidean QCD (and other gauge theories) can be computed exactly in a particular finite- 
volume scaling region. The origin of these developments dates back to work in the 1980's on QCD in 
a finite volume (see e.g. ref. Q), but the main breakthrough came with two very influential papers 
in 1992 by Leutwyler and Smilga [Q], and by Shuryak and Verbaarschot |J. In the latter paper an 
intriguing relation to random matrix theory was pointed out for the first time, and this led quickly 
to a series of theoretical developments that clarified the connection between Dirac eigenvalue spectra 
in gauge theories and random matrix theory ||, |6| (see also the reviews of ref. [JjJ). Here we shall 
discuss some very recent developments. 

The central object of study is the spectral density p(A; mi, . . . ,rriN f ) = J2n(^(^ ~ ^n))u, where the 
average is taken over all gluon configurations with fixed topological charge u, and where the Dirac 
eigenvalues A n are solutions toI)(p n = \ n 4> n . In a finite volume V it is convenient to introduce instead 
a rescaled, double-microscopic, spectral density ||] 

ps((^i,---,m f ) = > v - 00 > C 1 ) 

which computes the local density near A ~ 0, on a scale set by the chiral condensate X, as computed 
in the massless theory. In the infinite- volume limit X will be proportional to p(0), but the microscopic 
spectral density will generically vanish at £ = 0, in accordance with the fact that in any finite volume 
V there is no spontaneous chiral symmetry breaking. 

Consider now a finite- volume range « L « l/m n , where L ~ V 1 ^. The essential 

observation of Leutwyler and Smilga was that in a sector of fixed topological charge v this actually 
defines a finite-size scaling region. What this means becomes clear when one considers the finite- 
volume partition function (really the generating function for the chiral condensate). In the above 
limit it equals 

zi Nf) Qui,..., fi Nf ) = J dU(det U) v exp [VE Re Tr[MU]] (2) 

where the integral is taken over U(Nt), and A4 is the quark mass matrix, which we take to be 
diagonal in the masses m^. This generating function depends only on one very particular combination, 
Hi = niiVT,, and is in this sense a scaling function. The only needed ingredient is the existence of a non- 
vanishing chiral condensate X. Although the effective Lagrangian is consistent with the Gell-Mann- 
Oakes-Renner relation, the partition function does not depend on other dimensionful parameters, such 
as fn. Of course, the above representation becomes exact only in the limit, but this is precisely what 
we mean by having an exact finite-size scaling function: the universal result can be recovered to any 
required accuracy by tuning V. Corrections are suppressed by powers of 1/V; we will comment on 
the nature of such corrections below. 

The finite-volume partition function (||) has been evaluated exactly for an arbitrary mass matrix M. 
and for any Nj and v H. It turns out that complete spectral information about the Dirac operator 
in the double-microscopic limit can be obtained from this effective partition function alone H [To| . 
Curiously, the simplest starting point is the random matrix theory formulation of the finite-volume 
partition function @ ||: 



zl Nf \mi,...,m Nf ) = dW Y[det{M + m f ) exp 



—trV(M 2 ] 



2 



(3) 



1 



where 




Here W is a rectangular complex matrix of size N x (iV+|f|). In the large- iV limit the space-time 
volume V of QCD is identified with 2N. The potential V(M 2 ) can be parametrized in a general 
way by V(M 2 ) = J2(9k/k)M 2k . It was proven in refs. |], |ll[ that all double-microscopic spectral 
correlators (including, as the most simple case, the double-microscopic spectral density itself) in fact 
are universal, i.e. independent of the choice of V(M 2 ) up to a rescaling of the local macroscopic 
spectral density, here p(0). The class of potentials that fall into this specific universality class is huge, 
its boundary given by potentials for which one has p(0) = 0, but p( 2n \0) for some integer n |E 



To begin, the random matrix theory computations were always performed in the specific case of a 
Gaussian potential. Rewriting that particular random matrix theory partition it was shown already in 
refs. H [13|] that in the double-microscopic limit in which £ = AA r 27rp(0) and m = niiN2irp(0) are kept 
fixed as N —* oo this Gaussian random matrix partition function equals the Leutwyler-Smilga chiral 
Lagrangian (|2|) provided one makes the identification £ = 2-7rp(0). This exact equivalence between 
the random matrix partition function and the QCD effective Lagrangian in this regime, eq. (||), is 
crucial for the understanding of why random matrix theory can be used to compute Dirac operator 
spectra. It is therefore important to notice that also this equivalence holds universally, independent 
of the choice of the random matrix potential (up to the restrictions specified in refs. |6], |l2|). We shall 
now prove this. 



2 Universality of the partition functions 

There are at least two ways to demonstrate universality of the random matrix partition function, and, 
subsequently, the identity (up to an irrelevant /^-independent constant): 



(pi,...,p Nf ) = 2,i Nf \pi, . . . ,p Nf ) , (5) 



Here the l.h.s. is the finite- volume QCD partition function (||), and the r.h.s. is the random ma- 
trix partition function (||) evaluated in the double-microscopic limit. The most direct way is to 
use the orthogonal polynomial representation of the latter, i.e define monic orthogonal polynomials 
P n {X; mi, • • • , m N f ) so that (expressed in terms of A = A 2 , which are more convenient variables in the 
random matrix context): 

dX A^JJ(A + mj)e~ w(A) Pfc(A;mi,...,m A r / )P£(A;mi,...,m A r / ) = h k (m 1 , . . . , m Nf )5 k e • (6) 
/ 

The random matrix partition function is directly related to the normalization constants /*!&: 

AT-i JV-l 
zi Nf) = N\ Y[ h k (mi, . . .,m Nf ) = N\h (mi, . . .,m Nf ) N r k (mi, . . .,m Nf ) , (7) 

k=0 k=l 

where r k = h k /h k _i. It is now possible to use the universality proof of ref. ||, extended to the 
case of finite masses flllfl , to prove the universal relation (|5[) in the double-microscopic limit. This is 
most easily done by taking the logarithm of (0), and turning the resulting sum into an integral in the 
large- N limit. 
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A much simpler way to prove universality of (|5|) is to make use of an interesting relation between the 
orthogonal polynomials of eq. (||) and the random matrix partition functions [0]: 



id I 2 \ fit *t\N / \-v ^ U*l> • • • iPNf+l) , Q s 
^V(.-MjV J +l:A f l ! --- J MAT / J = ^(-J-J (MiVf+lJ Hv") • (°) 



Here C is an irrelevant /^-independent constant that just fixes the overall normalization of the poly- 
nomials. The l.h.s. of eq. (||) was proven to have a universal double-microscopic limit in refs. j|, 11] 



(and the solution is unambiguously defined also for negative eigenvalue entries). The universality 
proof now proceeds recursively, or by induction. For Nf = (the quenched case) universality holds 
trivially. Universality of the random matrix partition functions for Nf = 1, and higher, then follows 
recursively, using eq. (|8j). Of course, this proves only universality of the result from the random matrix 
side, but not the identity (|5). Fortunately, the identity has been established for the special case of 
Gaussian potentials in refs. |3|, [l~3[] . The simple proof given here therefore extends this identity to the 
full universality class. 

3 Spectral correlators from partition functions 

The very useful connection between the orthogonal polynomials and random matrix partition functions 
(||) is actually only one in a series of such relations 0, 14 1. Taken separately, these relations provide 



very convenient and compact expressions for all the relevant objects that enter in the random matrix 
computations. But when used in connection with the identity (||) they provide a much more intriguing 
series of relations - relations that now only refer to the finite-volume QCD partition functions. 

Most important in this context is the corresponding partition function representation of the kernel in 
random matrix theory: 

N-l 

K N (X,X'; mi ,...,m Nf ) = C "T C^ 2 W 2 )) (AA > r +§ TJ ^ (A 2 + m 2 )(A « + m 2 } £ P,(A 2 )P,(A' 2 ) . 

/ i=0 

(9) 

As is well known, from this kernel one can derive all spectral correlation functions in the limit TV — » oo: 
p(Xi, . . . ,X n ;m 1 , . . . ,m N ) = det K(X a , X b ; mi, . . . , rriN f ) ■ (10) 

' a,b ' 



We can now make use of the following convenient representation of the kernel: 



K N (z,z';mi, ■ ■ -,m Nf ) 



-N {v{z 2 ) + V{z ,2 )) + l Nf 



e 2 



{zz') u+ 2 n V/ J(z 2 + m 2 Mz> 2 + m 2 fi 



N f r ooN-lf N, 

x riK) / n UhKih - ^ 2 )(a, - ^ ,2 )ii( A * + 

/ i=l \ / 



Z v (mi, . . . ,mN f ) 



m})e- NV ^ 



detijXy 1 



(11) 



Except for the fact that the last eigenvalue integral runs up to — 1 only, the last factor is simply 
yet another partition function, now with two additional quark species of imaginary mass! Thus, up 
to corrections of order 1/N, we have in the large- A limit: 



N f 

KZ f '"\z,z!;m 1 ,...,m N/ ) = e~% ^ z ^ +v ^ (-iy Vz^J{ ^(z 2 + m})(z> 2 + m)\ 
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Z u (mi, . . .,m Nf ,iz,iz ) 

x — jwr ( , — ' (12) 

Z u (mi, . . .,m Nf ) 

We are now ready to take the double-microscopic limit in which £ = zN2irp(0) and pi = miN2irp(0) 
are kept fixed as iV^oo. In this limit the prefactor exp [-(N/2)(V(z 2 ) + V{z' 2 ))} becomes replaced by 
unity. By identifying £ = 2-7r/)(0), and using the universal relation (||) we finally arrive at the following 
master formula ||: 

N f (N f +2), 

4^(C,CW..,/^) = CtVTYl + + - , (i3) 

/ Zv {pi, ■ ■ -,/J-Nf) 

From this one single formula all double-microscopic spectral correlators can be computed directly from 
QCD chiral Lagrangians in the appropriate scaling regime. In particular, for the spectral density itself 
we find 

/MC/M = C 2 |Cin(C 2 + ^) Z ^r'''' ,m/,<,<) - (14) 

/ Z u (p, t , . . . , p, Nf ) 

The overall proportionality factor C2 can be determined by using the matching condition 



lim p s (C;m,...,fi N .) = 1/tt , (15) 

Q— >CJO 



which fixes C 2 = (-iy+t N f/ 2 l 



The higher /c-point double-microscopic spectral correlation functions are conveniently evaluated using 
the double-microscopic limit of the general relation ([To|). Curiously, it is also possible to relate these 
higher fc-point functions to finite-volume QCD partition functions with 2k additional quark species 



Ps ((i,...,(k;»u...^N f ) = (CiII(C? + M/) I nici-c 



2|2 



zj Nf+2k \p 1 , . . .,HN f ;{i&}, ■ ■ • ,{»Cfc}) 

Zy (/Ul, . . .,PN f ) 



(16) 



Each additional imaginary quark mass iQ is thus doubly degenerate. The overall proportionality 
constant can again be fixed by a matching condition. For k = 1 the relation ( |l~6| ) simply coincides 
with the previous expression for the double-microscopic spectral density. But already for k = 2 (and all 
higher values of k) the expressions are completely different, relating as they do the spectral correlators 
to finite-volume QCD partition functions with different numbers of flavors. It is quite amazing that 
the finite-volume QCD partition function (|2|) has all this structure, which takes on such a simple form 
in random matrix language, encoded in it. In fact, by combining eqs. (|l0|) and ( |l6| ) one obtains 
an infinite sequence of consistency conditions for QCD partitions. The relations become particularly 
transparent if we first take the additional fermion masses to physical values by replacing Q — ► — iQ 
(inspection of the explicit solution of ref . || shows immediately that this can be done unambiguously) . 
We then find the following infinite sequence of consistency conditions [0]: 



det 

l<a,fe<fc 



VCZ II vM-C 2 )(^-C^f /+2) (mi , • • • ,m f , (a, C&) 
/=i 
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k ( Nf \ k gr( N f+ 2k )( it \ if \\ 

n c, n w - c 2 ) im - <? i 2 - — • (17) 

where is some overall /^-independent normalization constant. Precisely these relations encode in 
the finite-volume QCD partition function the fact that in the random matrix picture the kernel (||) 
generates all spectral correlation functions through the relation (|i~0|). 



4 Direct computations from chiral Lagrangians 



We have learned that the massless spectral sum rules § do not provide the proper starting point for 
computing the microscopic spectral density. It is the double-microscopic limit 



15, 11 



] that is 

needed. This was in fact clear already from the first demonstration of the equivalence of the random 
matrix theory partition function and the finite- volume QCD partition function ||, 13]. Instead of 
the massless spectral sum rules, one should focus on the "massive spectral sum rules" || |l^] because 
these contain the analytical structure that allows one to unravel the spectral correlators from the 
finite-volume partition function. This fact becomes very clear when one considers the most simple 
example, the massive spectral sum rule corresponding to quenched QCD. Defining 



G(n) = 2fj, / dX 
Jo 

this can be written as a Stieltjes transform: 



PsW 

A 2 + fJL 2 



°° dt Ps(t)/Vt 
o t + y 



G(y/y) = F{y) . 



The inverse of this is given by the discontinuity: 

^ = ilnn[F(-t 



F(-t + ie) 



(18) 



(19) 



(20) 



So if one could compute the l.h.s. of ( |l8|) directly from a finite-volume partition function, one would 
have achieved a derivation of the spectral density without having at any intermediate stage to go 
through the random matrix theory framework at all. The trouble is that the massive spectral sum 
rule (|l^) refers to a "quenched quark" of (rescaled) mass [i. If it were a dynamical quark, one could 
compute the function G{p) straightforwardly from [0 

G{fj) = — lnZ^/i)- - , (21) 



where the last term subtracts the contribution from the zero modes. 



The needed trick, recently 
discovered by Osborn, Toublan and Verbaarschot [10], is to compute the r.h.s. of eq. (21) in a 
finite- volume field theory that contains yet another "quark" , now of opposite statistics and of initially 
different mass (so that, after taking the degenerate mass limit, the two determinants cancel in the 
partition function itself). The result is (for v = 0) ]lO| ]: 



(22) 



a result that was first derived the other way around, from the random matrix theory result, by 
Verbaarschot | |18| ]. Substituting this into eq. (p0|), one finds, straight from the finite- volume partition 
function, 

1, 



PS(A) = -|A| j Q (\y + Ji(\ 2 ) 



(23) 
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which of course agrees with the result obtrained from random matrix theory. Not only could those 
authors compute the function G(n) this way, they also managed to rewrite the general expression (|20| ) 
in precisely the form (|TJ) using the technique of partially quenched chiral perturbation theory |19|] , 
here based on the super Lie group U(Nf + 1|1). Their result naturally generalizes to higher fc-point 



spectral correlation functions, now given in the form (16). The relevant super Lie group will here be 
U(Nf + k\k). 



5 Finite-volume corrections 

We have seen that the microscopic spectral correlators have a natural interpretation as finite-size 
scaling functions. This makes them ideally suited for lattice gauge theory studies, and in fact there 
have now been Monte Carlo tests of the spectral densities of QCD in both (3+1) dimensions |2l|] 



and (2+1) dimensions [22]. (There have also been interesting studies of the applicability of random 
matrix techniques beyond the microscopic limit, in the "bulk" pq]). One obvious question in that 
connection concerns finite-size corrections. In actual computations the volume V is often far from being 
asymptotically large, and one could ask whether it is also possible to analytically calculate subleading 
corrections. For example, for the double-microscopic spectral density itself one could envisage an 
expansion of the kind 



P { P(CiV1, ■ ■ -,V>N f ) = Ps"'(.C,Pl, ■ ■ -,VN f ) 



(oo) . 



i + T?/(C;m, • • -,m f ) + • 



(24) 



with A some dimensionful constant, and / a correction-to-scaling function. One could hope that such 
corrections could be computed analytically using the random matrix theory formulation. In fact, if we 
go back to the derivation of eq. (12) from (|Tl|), we could try to keep the subleading corrections that 
come from ignoring the difference between N and N — 1 in (|ll]). However, on top of these we can also 
get subleading contributions from the potential V(X 2 ), even in the microscopic limit. We therefore 
conclude that such subleading 1/N corrections in the random matrix picture will be non-universal, 
and hence cannot be expected to be related to the Dirac eigenvalue spectrum. This is completely in 
accord with the field theory picture, in which subleading terms in l/V will involve non-static modes 
of the pseudo-Goldstone bosons. The kinetic term ^f^Trld^Ud^U] in the effective Lagrangian can 
therefore not be neglected. A new dimensionful scale (namely f n ) has entered, and the Dirac spectrum 
will cease to be a scaling function related to just V and E. Of course, one could try to systematically 
analyze l/V corrections in this very precise framework of the effective Lagrangian. The most natural 
starting point will unfortunately not be the usual expansion around the kinetic term, but rather a 
low-momentum expansion around the mass term TrLVli/]. 



6 Flavor dependence 

We finally address the question of the flavor dependence of all these results. The number of flavors 
enters in a very simple way in both the field theory and random matrix picture. In the former it 
determines the coset integration for the effective Lagrangian, while in the latter it enters only through 
the strength of the determinant in the expression @. Both lead to an extremely mild dependence 
on the number of flavors. This is because we throughout normalize the chiral condensate to one 
flavor-independent number S, - a convenient normalization because it puts the different theories with 
different flavor content on the same common scale. Clearly the whole framework collapses as the 
number Nf exceeds the value NJ above which QCD no longer supports spontaneous chiral symmetry. 
If we allow ourselves to treat this upper number of flavors N* as a free and tunable parameter, then 



6 



a normalizable condensate S will simply cease to exist precisely at Nf = Nf. In the random matrix 
theory context this may correspond to hitting the boundary where p(0) — * 12]. Beyond this point all 
results discussed here will no longer be valid. Just at the point where the condensate disappears one 
can define critical exponents that count the rate at which p(X) vanishes as A — ► 24]; however they 



will here not correspond to a physical phase transition (appearing instead as one tunes Nf continuously 
to Nf). Beyond that point the spectral density of the Dirac operator will develop a gap around A = 0, 
and there is no obvious way to compute it from the finite-volume partition function (mesons with 
quantum numbers of the pseudo-Goldstone bosons will still be the lightest excitations pq| , but there 
will be no analogue of the chiral Lagrangian). 

There is however a limit of large Nf in which the results here will continue to be valid. This is what 
is called the "topological 1/N expansion", where Nf — > oo and N c — * oo, with the ratio ij = Nf/N c 
kept fixed. The pertinent random matrix theory ensembles (and chiral Lagrangians) continue to be 
the same, and if r\ is chosen small enough (as in QCD), the theory will still undergo spontaneous 
chiral symmetry breaking. However, because Nf — > oo, there is now non-trivial "dynamics" even in 
the very simple chiral Lagrangian (||). In fact, when all masses are chosen equal the theory becomes 
identical to large-iV c lattice QCD in (1+1) dimensions will so-called Wilson action. It is known that 



this theory undergoes a 3rd order phase transition [26], and the QCD partition function (g) therefore 



undergoes precisely such a phase transition in the above limit. Because of the connection between 
partition functions and the microscopic spectral correlators, such a phase transition is expected to 
surface also in the Dirac spectrum, once the above limit is taken. Because a new limit (Nf — > oo and 
N c — ► oo, with r\ = Nf/N c fixed) is taken on top of the usual double-microscopic limit, we need to 
redefine our scaling variables. Let k = Nf/fx, and keep this variable fixed as Nf — > oo. If we define 
the free energy by T = — [In Z( N f']/Nf, then from the Gross- Witten analysis we know that 

T = 1/(4k 2 ) , k > 1 

T = 1//c + (1/2)1ii(as)-3/4 , k<1. (25) 

Taking derivatives, one indeed verifies that only the 3rd derivative of T is discontinuous at k = 1. It 
is amusing that this 3rd order phase transition of the effective QCD partition function has a simple 
interpretation as occurring in a two-dimensional "world volume" taken to be the 2-dimensional space 



of the Nf x Nf unitary matrices U [27]. One can plot the complex zeros of the partition function, 
and see that they as expected precisely pinch the real K-axis at K = 1. From this one can very 
accurately compute the correlation length critical index v, which turns out to agree with an analytical 
result obtained from the so-called double-scaling limit of models of 2-d quantum gravity pq] . More 
importantly, the critical scaling corresponding to this critical index sets in at very, very low values 



of Nf - actually already at Nf = 2 [27|. Traces of this analog of the Gross- Witten phase transition 



therefore persist even in real (3+l)-dimensional QCD, in the microscopic scaling regime. 

Acknowledgements: I thank G. Akemann, U. Heller, U. Magnea, S. Nishigaki, and T. Wettig for 
stimulating discussions in connection with our joint work. 
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